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What is claimed is 

5 A method for calculating position and time of a GPS 

receiver ccW>rising: ' ^ jLidfr^^ 

providing pseudo ranges that^esti-ma-te- the range of the 
GPS receiver tfe a plurality of GPS satellites; 

providing an estimate of an absolute time of reception 
0 of a plurality of Nsatellite signals; 

providing an estimate of a position of the GPS 
receiver ; 

providing satellite ephemeris data; 

computing absolute position and absolute time using 
said pseudoranges by updating said estimate of an absolute 
time and the estimate of position of the GPS receiver. 



2 . The method of claim 1 wherein said pseudoreanges are 
sub-millisecond pseudoranges. 

0 

3. The method of claim 1, wherein a first instance of said 
estimate of absolute time is in error by more than 10 
milliseconds . 

5 4. The method of claim 1, wherein said estimate of absolute 
time is provided by a clock that is not linked to a GPS 
reference time. 

5. The method of claim 1, wherein said estimates of 
0 position and absolute time may be arbitrary guesses. 



6. The method of claim 1, further comprising: 

forming a-priori pseudorange-residuals ; and 
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forming a mathematical model relating said a-priori 
pseudorange-residuals to said updates of said position and 
absolute time estimates; and 

computing updates of the position and absolute time 
5 that fit the mathematical model. 
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7. The method of claim 6, wherein said a-priori range- 
residuals are a difference between expected ranges, from 
said satellites to said a-priori position estimate, and the 

1 0 pseudoranges . 

8. The method of claim 6, wherein said expected ranges are 
computed at a time given by said a priori time estimate. 

15 9. The method of claim 6, wherein said mathematical model 
is a linearization of a Taylor series of a non-linear 
mathematical model . 

10. The method of claim 9, wherein said linearization is of 
2 0 the form: 

w, = / dx, dp, I dy, dp i I dz, 9p, / dt c , dp, I dt s ]* [x, z, t c , t s J 



where : 

25 Ui is the a-priori range-residual, for one satellite; 
p i is the pseudorange for a satellite i; 
x, y and z are three coordinates of position updates; 
t c is a common mode error update; 

t s is an absolute time of reception update; and 
3 0 d denotes the partial derivative. 

11. The method of claim 4, wherein one or more of the 
updates is assumed known, so that the remaining updates may 
be computed. 
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12. The method of claim 10, wherein one or more of the 
updates x, y, z, or t c is assumed known, and set equal to 
an assumed value in the mathematical model, so that the 
remaining updates may be computed. 

5 

13. The method of claim 10, wherein one or more of the 
updates x, y, z, or t c is assumed known, and added to the 
model as a pseudo-measurement so that the remaining updates 
may be computed. 

0 

14. A process as in claim 6, wherein other measurements or 
constraints are used in the mathematical model . 

15. A process as in claim 6, wherein said a-priori 

5 position is obtained from a location of a radio tower used 
to communicate with a mobile device containing said GPS 
receiver . 

16. A process as in claim 6, wherein said absolute time 
0 estimate is obtained from a real time clock in a server, 

said server being located remotely from said GPS receiver. 



17/ A metis^d for calculat ing a bsolute ti me of a GPS 
receiver composing: 

providing p^eudoranges that e stimate the range of the 
GPS receiver to a^plurality of GPS satellites; 

providing an estimate of position of the GPS receiver; 

and 

computing absolute t^jjie using the pseudoranges and the 
3 position estimate. 

18. The method of claim 1 wherein said pseudoranges are 
sub-millisecond pseudoranges. 
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19. The method of claim 17, further comprising: 

forming a-priori pseudorange-residuals ; and 
forming a mathematical model relating said a-priori 
5 pseudorange-residuals to updates of said estimate of 

position and absolute time; and 

computing updates of the position and absolute time 

that fit the mathematical model. 
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10 20. The method of claim 19, wherein said a-priori range- 
residuals are a difference between expected ranges, from 
said satellites to said a-priori position estimate, and the 
pseudoranges . 



15 21. The method of claim 19, wherein said expected ranges 
^ are computed at a time given by said a priori time 

estimate . 



a 22. The method of claim 19, wherein said mathematical model 

20 is a linearization of a Taylor series of a non-linear 
mathematical model . 



23. The method of claim 22, wherein said linearization is 
of the form: 

25 

= [3p, / dx, dp i I dy, 3p, / 3z, 9p, / dt c , dp, I dt s ]* [x, z, t c , t s J 

where : 

u i is the a-priori range-residual, for one satellite; 
30 p ± is the pseudorange for a satellite i; 

y and z are three coordinates of position updates; 
t c is a common mode error update; 

t s is an absolute time of reception update; and 
d denotes the partial derivative. 
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24. The method of claim 23, wherein one or more of the 
updates x, y, z, or t c is assumed known, and set equal to 
an assumed value in the mathematical model, so that the 
remaining updates may be computed. 

25. The method of claim 24, wherein one or more of the 
updates x, y, z, or t c is assumed known, and added to the 
model as a pseudo-measurement so that the remaining updates 
may be computed. 

26. A process as in claim 24, wherein other measurements 
or constraints are used in the mathematical model . 

27. A process as in claim 24, wherein said a-priori 
position is obtained from a location of a radio tower used 
to communicate with a mobile device containing said GPS 
receiver . 



A method of calculating a GPS position for a GPS 
receiver from partial pseudoranges that have ambiguity in a 
number of integer milliseconds, comprising: 

a) choosing an a-priori position of the GPS receiver; 

b) calculating integers conforming to said a-priori 
position; 

c) calculating a navigatio n sol ution; 

d) calculating a-posteriori residuals; and 

e) using a relative size of said a-posteriori 
residuals to determine if said calculated integers are 
correct; and 

f ) repeating steps c) , d) and e) using another a- 
priori position until residuals having a magnitude below a 
predefined threshold are computed. 
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29. A method as in claim 28, wherein said a-priori 
position is not within a specified distance of the actual 
receiver position. 

30. A method as in claim 28, wherein said a-priori 
position is greater than 100km from the actual receiver 
position . 

31. A method as in claim 28, wherein said a-priori 
position is greater than one integer millisecond from the 
actual receiver position. 

32. A method as in claim 28, wherein said a-priori position 
is an arbitrary guess. 



33i AXsystem for computing position and time for a GPS 



received comprising : 

a mol&ile device comprising a GPS receiver and a 
wireless transceiver ; 

a serveV being in wireless communication with said 
mobile device \ 

where said GPS \eceiver computes pseudorang es that estimate 
the range of the^G ps receiver to a plurality of GPS 



satellites and theV^ireless transceiver transmits said 
pseudoranges to the^server; 

where the server \computes an absolute position and 
absolute time for the Gl?S receiver using the pseudoranges 
and an estimate of positron and time. 

34. The system of claim 33 wherein said pseudoranges are 
sub-millisecond pseudoranges. 

35. The system of claim 33 wherein the position estimate 
is a position of a radio tower coupled to the server that 
receives signals from the wireless transceiver. 





